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ABSTRACT 

Context. Hard X-ray emission of coronal sources in solar flares has been observed and studied since its discovery in Yohkoh observa- 
tions. Several models have been proposed to explain the physical mechanisms causing this emission and the relations between those 
sources and simultaneously observed footpoint sources. 

Aims. We investigate and test one of the models (intermediate thin-thick target model) developed on the basis of Yohkoh observations. 
The model makes precise predictions on the shape of coronal and footpoint spectra and the relations between them, that can be tested 
with new instruments such as RHESSI. 

Methods. RHESSI observations of well observed events are studied in imaging and spectroscopy and compared to the predictions 
from the intermediate thin-thick target model. 

Results. The results indicate that such a simple model cannot account for the observed relations between the non-thermal spectra of 
coronal and footpoint sources. Including non-collisional energy loss of the electrons in the flare loop due to an electric field can solve 
most of the inconsistencies. 
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1. Introduction 

How well do we understand the physics behind solar flares? 
Solar flares release a large portion of their energy within sec- 
onds to minutes. A significant part of this energy goes into ac- 
celerated electrons and ions which then precipitate to the chro- 
mosphere along the field lines of a magnetic loop. The chromo- 
sphere acts as a thick target on the precipitating particles, lead- 
ing to characteristic hard X-ray (HXR) emission. The first ob- 
Ky servations of such footpoint sources were made by Hoyng et al. 
j_| ■ dl981| >. Masuda et al. (|1994| l found another HXR source n ear the 
& \ top of flare loops in Yohkoh observations. Feldman et al. d!994l > 
analyzed such coronal sources with Yohkoh BCS and SXT. Due 
to the normally low coronal densities one would expect a thin 
target emission from this coronal source. However, the obser- 
vations by Feldman et al. (1994) yield quite high column den- 
sities (10 20 crrT 2 < N < 7 ■ 10 20 crrT 2 ) in the coronal source that 
would act as a partly thick target on electrons below a certain 
energy. On the basis of the observations by Masuda et al. ( 1994 ) 
and Feldman et al. (11994K Wheatland and Melrose d 19951 ) de- 
veloped a simple model that has been used and investigated fur- 
ther (e.g. Metcalf & Alexander [19991 Fletcher & Martens [19981 
). The model consists of four basic elements; a particle acceler- 
ator on top of a magnetic loop, a coronal source visible in SXR 
and HXR, collision-less propagation of particles along the mag- 
netic loop and HXR-footpoints in the chromosphere. The coro- 
nal source acts as an intermediate thin-thick target on electrons 
depending on energy (thick target for lower energetic electrons, 
thin target on higher energies). We will therefore refer to this 
model as intermediate thin-thick target, or ITTT model. 

For large enough column depths or steep enough electron 
spectra, almost all electrons would be stopped in the loop and no 
footpoints could be observed. Such events were observed with 



Yohkoh and later with RHESSI. Veronig & Brown (2004) ana- 
lyzed flares with faint footpoint emission but a dense loop acting 
as a thick target on electrons of energies up to 60 keV. 

Due to its high spectral and spatial resolution, RHESSI 
(Lin et al. 2002) provides the means to perform imaging spec- 
troscopy on events with multiple sources. Battaglia & Benz 
(2006) showed that it is possible to study the non-thermal X- 
ray emission of coronal sources and footpoints and the relations 
between them in detail. 

Therefore, the ITTT model and its predictions can be tested 
using RHESSI data of well observed events. It is the simplest 
model that can explain the Yohkoh data. Can it also explain 
RHESSI data? 



2. Theoretical model 

The ITTT model features a dense coronal source into which a 
beam of electrons with a power-law energy distribution is in- 
jected. Some high energy electrons then leave the dense region, 
precipitating down to the chromosphere. The coronal region acts 
as a thick target on particles with energy lower than a critical 
energy E c and as thin target on electrons with energy larger than 
the critical energy. This results in a characteristic HXR spectrum 
as well as SXR emission due to collisional heating of the coro- 
nal region. The altered electron beam reaches the chromosphere, 
causing thick target emission in the footpoints of the magnetic 
loop. The predicted photon spectra are shown in Fig. Q] The 
coronal spectrum consists of two power-law components with 
a break at E c . The footpoint spectrum has a break at the same 
energy as the coronal spectrum and is a power-law at high ener- 
gies. 

The model predicts the following properties: 
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Table 1. Eventlist. The time indicates the peak time in the 50 
100 keV energy band. 



5 10" - 



Energy [keV] 

Fig. 1. Schematic photon spectra for electron distribution with 
power-law spectral index 6 = 5 and critical energy E c = 10 keV. 
The dashed line gives the spectrum of the coronal source with 
the spectral indices y°\ - y? - S-\ and - 5+1. The dotted 
line indicates the combined spectrum of the two footpoints. The 
footpoint spectrum is a power-law with spectral index y^ p =6—1 
at high energies. The solid line is the total of coronal source 
and footpoints (pure thick target). The dash-dotted line gives the 
footpoint spectrum from an electron population that has been 
altered in the loop (see Discussion). 



1 . The coronal and footpoint spectra intersect at the critical en- 
ergy (E c = E intersect ). 

2. There is a break at E c in the individual spectra of coronal 
source and footpoints, but not in the sum of the two spectra. 

3 . The spectral indices below and above the break in the coronal 
source have a difference of - y°i = 2. 

4. The spectral index of the coronal source below the break is 
equal to the spectral index of the footpoints above the break 

(rT = y{")- 

5. The difference between the spectral indices of the coronal 
source and the footpoints at high energies is yZ s - y^ = 2. 

6. Coronal source and footpoints have the same intensity at E c . 

7. SXR and HXR emission in the coronal source are spatially 
coincident. 



3. Observations 

3.1. Event selection and spectral analysis 

We analyzed the five events described in Battaglia & Benz 
(2006) in an interval of 3 RHESSI spin-periods (about 12 s) 
around the time of maximum HXR-emission. The events are 
listed in Table Q] 

Imaging spectroscopy applying the Clean algorithm for im- 
age reproduction was used to produce spectra of coronal sources 
and footpoints. Reasons for favoring Clean over Pixon are dis- 
cussed in Battaglia & Benz 2006 To compare the observations 
with the model, the footpoints where treated as one region and 
only a total spectrum over both footpoints has been computed. 
The coronal sources have been fitted with a thermal component 
and a non-thermal power-law, the footpoints only with a non- 
thermal component. Varying the starting fit parameters and the 



No. 


Date 


Time 


GOES class 


1 


4-Dec-2002 


22:47:02 


M2.7 


2 


24-Oct-2003 


02:48:42 


M7.7 


3 


l-Nov-2003 


22:33:14 


M3.3 


4 


13-Jul-2005 


14:14:30 


M5.1 


5 


30-M-2005 


06:32:06 


X1.3 



fit energy range provided a validation for the stability of the fit 
and an estimate of the range of spectral index values. The spectra 
and fits are shown in Fig. [2] 



3.2. Density and critical energy 

In the ITTT model, the critical electron energy for the transition 
between thick and thin target for a given column density No is 
given by equating the source size to half of the mean free path, 



E c = ^2KNq 



(1) 



with the constant K = 2ne A In A. The Coulomb logarithm InA is 
about 20 for the electron densities and temperatures in the pre- 
sented sample. The column depth the electron beam passes in the 
coronal source was computed using RHESSI and GOES data, 
following the treatment in Wheatland & Melrose ( 119951 ). For a 
source volume V, a total source diameter of V 1 ' 3 is assumed. An 
electron beam injected in the middle of the source would then 
travel half this distance, giving a column depth of: 



(2) 



where the electron density n e = y/EM/V using the observed 
emission measure EM and the volume V = A 3 ^ 2 . The area A 
was measured from the 50% contour of the coronal source in 
RHESSI Clean images at energy 10 - 12 keV, taking into account 
the Clean-Beam size and effects of the pixel size on the contour 
determination. 

The emission measure was computed from spectral fits to 
RHESSI full sun spectra and the spectrum of the coronal source 
only. As noted by Battaglia & Benz (2006), the thermal foot- 
point emission generally is very faint in RHESSI observations. 
Therefore, the thermal emission measured in full sun spectra is 
predominantly coronal emission. However, the spectral fittings 
deviate slightly, the temperature being lower and the emission 
measure higher for the imaging spectroscopy fit. We therefore 
use both as a range of possible column densities. 

Additionally, we computed the GOES emission measure 
and temperature for comparison. The range of all RHESSI and 
GOES emission measures then gives an estimate of the accuracy 
for the column densities and critical energies (Tabled. 

3.3. Position of coronal source at different energies 

The ITTT model predicts spatial coincidence of the coronal 
emission for all energies. To check for this, Clean images at dif- 
ferent energies where made and the centroid of the 50% emission 
of the coronal source was computed in order to get the position. 
At energies higher than 20 keV, the footpoint emission starts 
to dominate, making an accurate determination of the coronal 
source-position difficult. We determined the centroid positions 
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Fig. 2. Coronal (crosses) and footpoint (dots) spectra of the five analyzed events. The thin solid line indicates the thermal fit to the 
coronal source. The thick solid lines give the power-law fits to coronal source and footpoints. Their extension to low energies (not 
observed) is given as a dot-dashed line. The dotted line indicates the average critical energy. 
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Fig. 3. Images of all events at energies 10 - 12 keV. The 50% 
contours of the maximum of the footpoints are given (black). 
The white crosses indicate the position of the 50% centroid at 
energies 6-12 keV (thin) and 16-22 keV (thick) with error 
bars. 



for the energy range 6-12 keV (thermal) and 16-22 keV (par- 
tially non-thermal). 



4. Results 

The main numerical results are given in Table [2] We find col- 
umn densities between 2.1 • 10 19 < Nq < 1.2 ■ 10 20 crrT 2 . This 
yields a range of E c from 7.7 keV - 21 keV, significantly lower 
than the values calculated by Wheatland & Melrose for Yohkoh 
events (15 keV< E c < 40 keV). The comparison of the data 
with the predictions from the ITTT model is not always possi- 
ble due to the presence of the thermal component. Especially, 
the estimated critical energy is within the range dominated by 
the thermal component in all events. Comparing the data with 
the predictions from the ITTT model gives the following results, 
following the numbering in Sect. [2j 



1. 



The power-law fits of the non-thermal coronal and footpoint 
spectra intersect at energies E- mtersect - 7.6 - 23.1 keV. In 3 
events, this energy is within or just outside the estimate for 
the critical energy derived from eqns. ([D and ©. Therefore 
E c w E in1ersec1 . In two events (no. 2 and 4), E in1ersect is signif- 
icantly higher than the computed E c (E mtersecl > E c ). 
In the events with E c * Ei„, ersec/ the estimated break energies 
and intersection energies are in the energy range dominated 
by the thermal emission of the coronal source. A break in the 
coronal spectrum can therefore not be detected, but might 
still exist. In events nos. 2 and 4, where E- mtersect > E c , the 
intersection point of the power-law spectra is observed, but 
without a break, contradicting the prediction. 
In event 2, the intersection energy is well outside the thermal 
emission but y| s - + 2 in the coronal source, contrary to 
the prediction. 

In event 2, y c x s > y^ p , not equal as expected from the model. 

y" - y^ 2 P > 2 in three events and y" _ Y% < 2 in two events. 
In three events, E c and Ei ntersect of the non-thermal fits are 
within the thermal emission of the coronal source. In event 
2, Eintersect is observed with the coronal source being brighter 
than the footpoints at E c . 



4 
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Table 2. Summary of observed parameters to be compared with the theory. The power-law spectrum of the coronal source can only 
be observed at energies higher than the thermal emission, therefore no value for the part below the break could be determined. 



Event 


l 


2 


Y? 
7" 


3.5±0.7 


(6.3±0.1) 
6.3+0.1 


ri" f 


2.8±0.2 


2.5±0.1 


r 2 °-ri p 


0.7±0.7 


3.8+0.1 


No (cm- 2 ) 


(2.5-4.9)- 10 19 


(2.8-6.1)-10 


E c (keV) 


9.5-13.5 


10.3-15 


^intersect (keV) 


10.3 


23.1 



3 4 5 



4.2±0.4 5.6±0.2 5.9±1.0 

3.5+0.1 2.7±0.2 2.9±0.1 

0.7±0.4 2.9±0.3 3± 1.0 

(2.1-4.3)-10 19 (2.6-3.0)- 10" (7.1-12)10 

8.9-12.6 7.7-10.5 16.2-20.8 

7.6 20.2 20.9 



7. No significant positional differences between energies 6-12 
keV and 16 - 22 keV have been found (see Fig. [3J. Due to 
the increasing footpoint brightness at higher energies, an ac- 
curate determination of the coronal source centroid above 22 
keV is not possible. 

4. 1 . Energy input into the corona 

Wheatland & Melrose (1995 ) proposed that the collisional en- 
ergy necessary to produce the non-thermal bremsstrahlung emis- 
sion in the coronal source is sufficient to account for the ob- 
served heating of the coronal source. They referred to standard 
flare values (flux, duration, electron spectral index) to estimate 
the non-thermal energies. With the observations presented here 
the energy input into the individual coronal source can be deter- 
mined. 

We compute the change in thermal energy AE,i, in the coro- 
nal source between two times in the flare. The thermal energy 
is calculated according to E t / lerm = 3ksT y/EM ■ V with a filling 
factor of unity. 

The non-thermal energy input in the coronal source is com- 
puted from the non-thermal coronal spectrum. As cutoff en- 
ergy, we used the intersection point of the thermal and non- 
thermal emission. Saint-Hilaire & Benz (2005) found that the 
average relation between photon-turnover and electron-cutoff is 
^cutoff" ~ turnover ' 1-7- This gives an upper limit for the elec- 
tron cutoff energy and therefore a lower limit for the total energy 
input. 

Comparing the two lower limits, we find that the average 
non-thermal energy input is at least of the same order as the 
change in thermal energy i.e. AE„ „_,i, > AE,/, in agreement with 
the prediction of the ITTT model. For a lower electron cut-off 
energy, the non-thermal energy input could be up to an order of 
magnitude higher. 

5. Discussion 

The ITTT model as presented by Wheatland & Melrose ( 1995) 
predicts a critical electron energy E c marking the transition be- 
tween thin- and thick target emission. This energy is likewise the 
intersection energy Ej ntersect between the non-thermal spectra of 
the coronal and footpoint sources. How can it then be explained 
in the frame of the ITTT model that in events 2 and 4 the esti- 
mated E c is significantly lower than the observed Ej„i ersect l 

E c as calculated here may be a lower limit. Including col- 
lisional deflections of the electrons within the coronal target 
rather than just energy-loss, the electron paths would become 
longer and the critical energy for thin target emission higher. 
This would increase E c and could bring it up towards the ob- 
served Ei n t e rsect- In that case however, one would expect to ob- 



serve a break at Ei ntersect . This is not the case. Further, the differ- 
ence yl s - y 2 f p is higher than 2 in both of these events, making 
this explanation even less likely. 

An attractive possibility to interpret both, the absence of 
a break and the larger difference in spectral index is non- 
collisional energy loss of the electrons in the loop due to the 
electric field caused by the return current. This would result in a 
harder and fainter footpoint spectrum, shifting E mtersect to ener- 
gies higher than E c . The expected footpoint spectrum in the case 
of electrons loosing 15 keV of energy in the loop is indicated in 

The remaining discrepancy then is the difference in spectral 
index yl s - yi being smaller than two in events 1 and 3. Such 
a difference could result where the targets cannot be clearly sep- 
arated into thin and thick, for example if the sources are close 
as is the case in those events. Part of the emission selected for 
the footpoints could originate from the coronal source, giving a 
mixed thin-thick target spectrum for the footpoints. More obser- 
vations are necessary to confirm this explanation. 

An estimate of the non-thermal energy input into the coro- 
nal source compared to the change of the thermal energy con- 
tent in the coronal source shows that collisions as manifest in 
hard X-ray emission may account for the observed heating of 
the coronal source. Depending on the electron cut-off energy, 
the amount of energy deposited in the target might even exceed 
the observed thermal energy. However, we did not consider cool- 
ing of the coronal source due to thermal conduction or radiation. 
Such effects could lead to a thermal energy input larger than ob- 
served. 

6. Conclusions 

We showed that a simple thin-thick target model as the one pro- 
posed by Wheatland & Melrose ( 1995) cannot explain all obser- 
vations made with RHESSI. Modifications of the model are nec- 
essary, the simplest being the consideration of the electric field 
due to the return current. In the absence of a balancing ion flux, 
such a return current is predicted from basic physics in view of 
footpoint sources. 
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